Exact density functional for hard rod mixtures derived from Markov chain approach 



(N 

o 

(N 
S3 

a 



in 



i 

S3 

o 
o 



> 

O 

o 



X 



Benaoumeur Bakhti, Stephan Schott, and Philipp MaastQ 

Fachbereich Physik, Universitat Osnabriick, Barbarastrafle 7, 49076 Osnabruck, Germany 

(Dated: January 19, 2012) 

Using a Markov chain approach we rederive the exact density functional for hard rod mixtures 
on a one-dimensional lattice, which forms the basis of the lattice fundamental measure theory. The 
transition probability in the Markov chain depends on a set of occupation numbers, which reflects 
the property of a zero-dimensional cavity to hold at most one particle. For given mean occupation 
numbers (density profile), an exact expression for the equilibrium distribution of microstates is 
obtained, that means an expression for the unique external potential that generates the density 
profile in equilibrium. By considering the rod ends to fall onto lattice sites, the mixture is always 
additive. 
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The extension of density functional theory from con- 
tinuum to lattice fluids [l| has proven to be useful for 
treating problems like ordering transitions [InSL proper- 
ties of interfaces separating different phases , phase 
separation in mixtures Q , or polymer adsorption at solid- 
liquid interfaces Q. Time-dependent density functional 
theory [H furthermore allows one to describe the kinet- 
ics of lattice fluids [loj . as emerging in phase ordering 
phenomena relaxation process es Jrj . and particle 
transport in driven lattice gases |13l4l5| . 

In 2002 Lafuente and Cuesta extended Rosenfeld's fun- 
damental measure theory to lattice models based on a 
derivation of an exact density functional for hard rod 
mixtures in one dimension [Tq . [l7j . This derivation was 
carried out following a procedure developed by Vander- 
lick et al. [l8[ for continuum fluids. Since the excess free 
energy part of the functional could be expressed in terms 
of differences between parts that agree in their functional 
form with the excess free energy functional of a zero- 
dimensional cavity, approximate functionals in higher di- 
mensions were obtained by dimensional expansion of the 
corresponding difference operator. By construction these 
fundamental measure functionals have the property to 
become exact under dimensional reduction and their im- 
pressive power was first shown by determining phase di- 
agrams of hard squares [13, EH and hard cube mixtures 
[lfjl ITU |20] ] with good quality. The fundamental measure 
functionals moreover allow one to apply the method of 
dimensional crossover and the merit of this was demon- 
strated by deriving functionals for lattice gases with near- 
est neighbor exclusion for different lattice types (square, 
triangular, face- and body-centered cubic) from the func- 
tional for cubes in (d+ 1) dimensions |21|. The structure 
of the corresponding results led to a suggestion how to 
construct fundamental measure functionals for hard core 
lattice gases for any type of lattice, shape of the particles, 
and arbitrary dimension [22| . 

In this report we rederive the exact density functional 
for hard rod mixtures in one dimension, that means the 
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starting point of the fundamental measure theory for 
hard core lattice gases, by applying the Markov chain ap- 
proach developed by Buschle et al. [23[ . This approach is 
conceptually different from the procedure of Vanderlick 
et al. [18] and we believe that it is useful and impor- 
tant on the following reasons: (i) The derivation of the 
functional becomes surprisingly simple. Making use only 
of the constraints of mutual rod exclusions, the relevant 
transition probability in the Markov chain is determined 
almost without any calculation, (ii) The transition prob- 
ability is (conditionally) dependent on a spatial region, 
where at most one particle can be placed, i.e. that of 
a zero-dimensional cavity. In this respect it reflects a 
property which turned out to be decisive for the gen- 
eralized construction of fundamental measure function- 
als by Lafuente and Cuesta 22]. (iii) The simplicity of 
the derivation suggests that it can be extended to hard 
rod mixtures with additional (thermal) interactions, (iv) 
The derivation yields also an explicit expression for the 
probability distribution of microstates for a given density 
profile. This means that in the present case an explicit 
expression for the "Mermin potential" is obtained, i.e. 
the unique external potential that would generate the 
given density profile in thermal equilibrium. In addition 
to these points we show that it is not necessary to con- 
sider non-additive mixtures when mixed parities of rod 
lengths are present (i.e. rods with both even and odd 
lengths in units of the lattice spacing). 

The mixture is considered to consist of q types of hard 
rods with length l a , a = 1, . . . , q in the presence of an 
external potential. It is convenient (although not neces- 
sary) to order the lengths according to l\ > Z 2 > . . . > l q , 
where different types of rods could have the same lengths 
due to different coupling to the external potential. The 
rods are located on a one-dimensional lattice with L sites 
and we set the lattice spacing equal to one. The lattice is 
defined in such a way that the ends of the rods coincide 
with lattice sites and we introduce occupation numbers 
n?, j = 1, . . . , L, a — 1, . . . , q, to specify the microstate 
of the mixture. If the left end of a rod of type a is at site 
j, then nf = 1, else n" = (here and in the following 
Greek superscripts refer to the type and must not mixed 
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FIG. 1. Illustration of the set of occupation numbers affecting 
the occupation of site k. Any placement of the left end of a 
rod of type a at the sites j with k — l a + 1 < j < k means 
that site k is covered by a part of this rod. This implies (i) 
that if a left rod end is at site k, all occupation numbers in 
the set {n°}n = {nf\l < a < q,k - l a + 1 < j < k - 1} 
must be zero, and (ii) that in the set {nf }t, = {n"\l < a < 
q, k — l a + 1 < j < k} there can be at most one occupation 
number with value 1. 



up with exponents) . The mutual exclusion of hard rods 
implies the constraint n^rij = for j = k, . . . , k + l a — 1 
(and k — j, . . . , j + lp — 1) 24]. In a grand-canonical 
description the chemical potentials \i a specify the mean 
numbers of rods of type a. 

To set up the Markov chain approach following [23[ 
it is useful to introduce the multicomponent state vari- 
ables fij — (nj, . . . ,n|) that can assume (q + 1) states 
e ,...,e g , where eo refers to an empty site, i.e. e"o = 
(0, . . . , 0), while e a , a = 1, . . . , q, refer to a site occupied 
by rods of type a, i.e. e a = (0, . . . , 1, . . . , 0) with the 1 
at the (a + l)th entry. The probability x(^ii • • • i ^l) of 
microstatcs can be decomposed as 



hard rod constraints, which reflects the corresponding 
property of a zero-dimensional cavity. In fact this set 
corresponds exactly to the zero-dimensional cavity for a 
mixture introduced in [TtJ as a collection of sets for each 
rod type. The property to have at most one occupation 
variable n" = 1 in {n"} k can be utilized to determine 
the conditional probabilities by simple probabilistic con- 
siderations. First let us write for a = 0, . . . , q 



u* - ir P\ \ ProbK = e a ,{n.} fc _i) 

ip[n k = e a \{n^} k -t) = J — (3) 

Prob({nj}fc_i) 

where Prob(.) denote joint probabilities. If a ^ 0, 
then all rij in the set {n^}/c-i must be zero. This im- 
plies Prob(n fc = e Q ,{n^} fc _i) = Prob(n fc = e a ,{n^ = 
0}fe-i) = Pf!, where p£ = (n a ) is the mean occupation of 
site k ((...} denotes an average over the microstate dis- 
tribution x(ni, . . . , Ul))- Since with the same reasoning 
Prob({n^ = 1 , all other = 0}fc_i) = Pj, we further 
have 



q k-1 

Prob({n? = 0}*_i) + X; E P? = 1 ( 4 ) 

p=ij=k-lf,+i 

due to normalization. Accordingly, we obtain for a =/= 



ip(n k = e a \{n?}k-i) 
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where we used one of the weighted densities (weighted 
mean occupations) [26| 



x(ni,...,n L ) = Y[tp( 



nk\n k -i, 



,ni) 



(1) 



fc=l 



where ip(.\ . . .) denote the corresponding condi- 
tional probabilities. To keep the notation sim- 
ple, we have labeled the starting of the chain, i.e. 
^(ni)ip(n2\ni)^(h3\fi2,ni) ■ .., by the same symbol "ip" 
(meaning in particular that ip(ni) is not a conditional 
probability). By using the Boltzmann expression for 
the probability of microstates in the grand-canonical 
equilibrium ensemble, i.e. x ex P[ — a ( u ? ~ ^a)nf], 
where j3 — 1/k^T is the inverse thermal energy and uf 
the external potential, it can be proven [251 ] that the 
conditional probabilities satisfy the Markov property 



ip(n k \n k -i,. ■ ■ , hi) = ip(hk\{nf}k-i) 



(2) 



where {n"}fc_i = < a < q, k — l a + l < j <k — 1} 

denotes the set of occupation variables, which have an 
influence on the occupation of site k, see Fig. [T] 

In the set {nf} k = {nf\l < a < q, k - l a + 1 < j < k}, 
i.e. all occupation variables involved in Fig. [TJ there can 
be at most one occupation variable n™ = 1 due to the 
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a — 1 j — l — 7n 



appearing in the lattice fundamental measure theory |16| . 
If n k — eo there are two possibilities: Either one ele- 
ment in {rij} k -i is one, or all elements are zero. In 
the first case, h k must be equal to eo, implying that 
the corresponding conditional probability is one. In the 
second case we need Prob(n k — eo,{f^ = 0}k-i) = 
Prob({n^ = 0}k) in Eq. ([3]), which by utilizing normaliza- 
tion as in Eq. Q (now with inclusion of site fc), is given by 



J - 



Prob({- 
In summary, 



ip(n k = e |{nf } fc _i) 



0} fc ) = i-ELiEi= fe -^ + i^ = i-5 r , 



(i) 

k 



(7) 



1 , one — 1 in {n?- } fc _i 



1 - S" (1) 

' ' all n? = in {nf } k -i 
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(0) 
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Combining Eqs. ([5]) and (J7]), we can write 

- 1 - Z-ip=l Z-,j=o n j 



1-5 



1 



j(0) 



n 



i-5 



(0) 



(8) 



where the distinction between the possible configurations 
in the set {rij }k is taken into account by the exponents. 

Inserting Eq. flSJ into Eqs. j2]) and (JTJ, the proba- 
bility distribution of microstates is given by the prod- 
uct of tp(hk\{rij}k-i) from Eq. flSJ) over all lattice sites, 
i.e. an explicit expression for x( n ) as function of the set 
n = {nf |1 < a < q, 1 < i < L} of occupation numbers is 
obtained (we define x( n ) — f° r au microstates n vio- 
lating the hard rod constraints). This means that, for a 
given density profile p = < a < q, 1 < < L}, 

the distribution of microstates is uniquely determined 
if we require it to satisfy the Markov property ([2]), i.e. 
x(n) = Xp(n). One could get the impression that this 
is more general than the uniqueness implied by the Mer- 
min theorem, which states that the prescription of p fixes 
the external potential u% — u^(p) in the sense that the 
Boltzmann distribution yields p in equilibrium in the 
presence of u£(p). However, since the Boltzmann dis- 
tributions satisfy the Markov property @, and Xp( n ) is 
unique, there is in fact no more generality, i.e. the mi- 
crostate distribution for given p satisfying the Markov 
property (j2J and the Boltzmann distribution generating 
p in equilibrium must be the same [27[. We can thus 
identify the "Mermin potential" U p (n) = J2 k u^(p)n^ 
by setting (3U p (n) oc — logXp(n), which, up to irrelevant 
constant contributions, yields (after some rearrangement 
of summations) 



k + l a 



<(p)=logp£-log(l-^ 0) )+ E l0 S 



1-5. 



(0)' 
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(9) 

Based on the Gibbs-Bogoliubov inequality the density 
functional in an external potential U(n) = J2k a u k n k IS 
defined as 



0(p) = ^ Xp (n) 



fc B riogx P (n) +U(n) -^TfJ, a N a 



a=l 



L q 



(10) 



k=l a=l 



where F(p) = k^,T^ n x P (n) log x P is the free energy 
functional. Inserting Xp( n ) one obtains 

L 

^F(p)=E{(l-4 1) )l°g(l-4 1) ) (11) 
(l-5< 0) )log(l-5f) + 5>*l°gPfc} 



fe=i 



Minimizing f2(p) with respect to the p" yields the density 
profile in equilibrium. 

Following Lafuente and Cuesta [17[, one can define an 
"ideal part" F id (p) by 



/3F id (p) = EI>feG°gp£-l) 



(12) 



k=l a=l 



This differs from the expression J^kiPk ^°EPk — (1 — 
J2a Pk ) l°g(l — E Q Pk )} f° r a non-interacting multi- 
component Fermionic lattice gas, but has the advantage 
to lead to a fundamental measure structure of the excess 
free energy part F cxc (p) — F(p) — F^(p). When using 



Eqs. ((nj, UT5J, and J2 a Pk = S k ~ S k this becomes 

L 

(3F cxc (p) = E{ [Si 1] + (1 - 5l 1} ) log(l - 5«) 
fc=i 

-[5f + (l-5f)log(l-5i°))]} (13) 

The terms in the square brackets have the same func- 
tional form as the excess free energy / cxc (ij) = r\ + (1 — 
r\) log(l — r\) of a zero-dimensional cavity with mean occu- 
pation r] [28[. Approximate fundamental measure func- 
tional in higher dimensions can be constructed by con- 
sidering the two terms in the square brackets as resulting 
from applying a one-dimensional difference operator and 
by generalizing this operator together with the weighted 
densities to higher dimensions (for details, see fl6l 

The excess free energy in Eq. (fTBf is equal to that 
found by Lafuente and Cuesta for an additive mixture. 
To recover their expressions, occupation numbers — 
0, 1 need to assigned to the rod centers, which amounts 
to a simple translation of the site indices, —¥ — 
n k+(i -e)/2' wnere e = if all l a are even and e = 1 if all 
l a are odd. 

Non-additive mixtures appear when considering a 
setup where the rod centers fall onto lattice sites and both 
even and odd l a are present, since in this case neighboring 
rods with even and odd l a have a minimum separation 
of half a lattice unit between their ends. For such non- 
additive mixtures one can construct the corresponding 
functional from that for additive mixtures [l7|. When 
the rod ends fall onto lattice sites, the mixtures are al- 
ways additive irrespective of having mixed parities of rod 
lengths. 
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